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ABSTRACT 



The purpose of this thesis is to investigate and establish 
Sturm-Liouvi 1 1 e properties for special eigenfunctions which 
are expressed in determinant form. In particular, a special 
case is presented where the elements of the determinant are 
Legendre polynomials. This type of determinant has a 
probability background dealing in birth and death processes. 
The method of analysis used in this thesis is a new approach 
to solving this specific example. This investigation involves 
systems of differential equations and Prufer's analysis in the 
phase plane. The following are new results obtained in 
addition to solving the special case mentioned above. 

Special determinants of hypergeometric functions also 
possess Sturm-Liouvil 1 e properties. As a special case, a 
different proof of Turan’s Inequality is provided. Finally, 
several theorems are presented for Sturm-Liouvi 1 1 e systems of 
differential equations with polynomial coefficients. 
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INTRODUCTION 



I . 



A . THE PROBLEM 

The main purpose of this thesis is to investigate and 
establish the essential Sturm-Liouvil le properties for special 
eigenfunctions which are expressed in determinant form. In 
particular, a special case is presented where the elements of 
the determinant are Legendre polynomials such that 



This type of determinant has a probability background 
dealing in birth and death processes and was studied by Dr . S. 
Karlin and Dr. J. L. McGregor. In their researches, Karlin 
and McGregor encountered determinants of classical polynomials 
as examples. In order to complete their results, they needed 
to know such properties as: number of zeros, interlacing 
zeros, and completeness. Their obvious starting point was to 
construct an ordinary differential equation such that (1.1) 
was the solution. However, this differential equation was not 
in a form that suggested how to answer any of these questions. 
[Ref. 1] 

The method of analysis used in this thesis is a new 
approach to solving this specific example. The investigation 
that follows provides a constructive answer, not only to this 




P}^{x) (n>k) k=0,l,2 

P„{x) Pjj^^(x) (-Kx<l) n=l,2,3 



# • • • 



/ • • • 



( 1 . 1 ) 
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example, but, to a whole class of such problems. The 
following are new results obtained in addition to solving 
Karlin's problem. 

Special determinants of other hypergeometric functions 
also possess the essential Sturm-Liouvi 1 1 e properties. 

As a special case of (1.1), a different proof of Turan's 
Inequality is provided in Chapter III. Turan's Inequality is 
of the form 



P^{x) 

P^*i(at) Pn*2(A:) 



< 0 



-Kx<l 



n=0, 1,2,... 



( 1 . 2 ) 



where the elements of the determinant are Legendre 
polynomials. 

Finally, several theorems are presented for 
Sturm-Liouvil le systems of differential equations with 
polynomial coefficients. In general , boundary value problems 
of systems are not well-posed. Until now, very little has 
been established concerning these systems with polynomial 
coefficients . 

B. STRATEGY FOR SOLVING THE PROBLEM 

The basic strategy developed here is to construct a second 
order differential equation such that (1.1) is a solution. 
Next, this differential equation is reduced to a special 
system of first order differential equations such that the 



2 



coefficients are polynomials. Then, the solution trajectories 
are analyzed in the corresponding phase plane by a method 
known as the "Priifer analysis". From this analysis, several 
theorems of Sturm-Liouvi 1 le type are developed for general 
systems of differential equations with polynomial 
coefficients. Finally, completeness and orthogonality are 
established with respect to a positive definite matrix weight 
function . 
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II. CONSTRUCTION OP A STURM-LIOUVILLE SYSTEM OF 
DIFFERENTIAL EQUATIONS 



A. DEVELOPMENT OP A NONSTANDARD STURM-LIOUVILLE DIFFERENTIAL 

EQUATION WITH SPECIAL DETERMINANT EIGENFUNCTIONS 

The objective of this section is to construct a second 
order linear differential equation, with appropriate boundary 
conditions, such that its solution is (1.1). For the general 
case, with arbitrary n and k, this method is somewhat lengthy 
and requires extensive manipulation of equations. However, 
considering a special case of k = 1, this method is not as 
complicated yet still illustrates all of the important steps. 
For completeness, the outline of the solution for the general 
case the same procedure and is provided in the Appendix. 

In the general case, the eigenfunction §jj]j(x) has the 
following form 



P,(x) 



(i2>ic) ic=0, 1, 2, . . . 

(-Kx<l) J3=l, 2, 3, . . . 



( 2 . 1 ) 



Setting k = 1 results in 






X %(3x2-l) 



PJx) p^^ix) 
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The strategy is to use recursion formulas and to develop 



three equations for §jj^(x), f'jjj(x), and f'jjj(x) as functions 
of PjCx) and P'jj(x). Then, from these equations, a system of 
equations in the matrix form AX = 0 is constructed where 
X # 0. Finally, the determinant of A is calculated and set 
equal to zero. This results in a differential equation having 
(2.1) as a solution. 

First, the equations for 5jjj(x), 
developed. Substituting the recursion formula 



Since k = 1, P|;(x) and P'|^(x) are replaced by x and 1, 
respectively. In addition, by computing the determinant and 
collecting like terms, we obtain the following equation for 




into (2.1) for P]j^j(x) and Pjj^ 2 (x) results in 



Pj, (X) xP* (X) - pI; (X) 










( 2 . 2 ) 
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Next, we differentiate both sides of (2.2). Thence, 






-X- 



‘Wi 



<^PiU)] . 



(2.3) 



Then, we substitute the Legendre differential equation 



d 

-ax 



(l-x^)P^(x) 



*=-7l(j3+l)P^(x) 



into (2.3) and collect like terms. Then, (2.3) becomes 



*i,,i(x) =(i3-l)xPa(x) + 



(2.4) 



Repeating the same procedure on (2.4) results in the 



following relation for f''uj(x) 



jp^(jc) + (72-l)XP^(x) 



(2.5) 



Now, a system of equations in matrix form is constructed 
from equations (2.2), (2.4), and (2.5). The system takes the 
form AX = 0 , i . e . , 
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AX= 



= 0 . ( 2 . 6 ) 







(n-l)x 


■ -1 • 


^n,lix) 


(23-l)X 


(n-1) (l-jc^) 


Pn(x) 


2(n+l) 




1-x^ 

2 


-xil-x^) 

n+1 


PLix) 



Since the trivial solution can not be a solution of (2.6), A 
is not invertible. Therefore, the determinant of A must equal 
zero. Hence, by setting the determinant of A equal to zero we 
obtain a second order differential equation where fjjj(x) is a 
solution . 



\A\ = (1-Jf2) (l+3x^)<,i(x) - 6x(l-JC*)*'.i(x) + 



(23-1) 



(3n+6)x^+23 






(2.7) 



(X) = 0 



By setting 1 = n-1, then (2.7) takes the form 



d-x") (l+3x»)<i(x) - 



3 (X+3)x^+A, 



-iK.i 



( 2 . 8 ) 



ix) 



Note that, for each Legendre polynomial solution P^(x) to 
the Legendre differential equation, there is a corresponding 
logarithmic solution, Q^(x). This logarithmic solution obeys 
the same recursion formulas as Pjj(x). Hence, by the same 
procedure, the second solution to (2.8) is 
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0„(x) Qo^iix) 



Therefore, the general solution is 





X 


(3x2-1) 




X %{3x2-1) 














Pn(x) 


Pa*i (X) 




Oa(x) Qa^y^(x) 



Finally, translated boundary conditions on (2.9) are 
imposed so the second logarithmic solution is eliminated. 

Note that (2.8) is not in the standard Sturm-Liouvi 1 le 
form. However, ijjj(x) still possesses all of the essential 
properties . 

B. REDUCTION TO A SYSTEM OF FIRST ORDER DIFFERENTIAL 
EQUATIONS 

The objective in this section is to reduce (2.8) to a 
system of differential equations of the special form 

(1-x^)®' = + (1+3 jc^)T -Kx<1 

( 2 . 10 ) 

(l-x^)T' = jg* + yT 



where a, fl, and y are polynomials in x, 5 = 5 jj 2 (x), and T is 
an arbitrary function in x. There are an infinite number of 
ways to transform (2.8) into a system of differential 
equations. However, by requiring polynomial coefficients we 
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remove the singularities at ± (-3)'*. Hence, all solutions are 
well behaved on the interval -Kx<l. 

The strategy is to convert (2.10) into a second order 
differential equation and then equate the polynomial 
coefficients of §* , §' and 4. 

The system (2.10) has the following equivalent form 



(1-x^) (3x*+l)®"- 



(a+Y+2jf) (3x*+l) +6jf(l-x*) - 



a' 0^*1) 1 * . 0 

1-X^ 1-X^ J 



First, note that the coefficients of 4* are the same. 
Next, equating the coefficients of 4' and simplifying results 
in 

y = -2x-a . (2.11) 



Equating the coefficients of 4 and substituting (2.11) for 
y leads to 



X 



3(X+3)x2+A-1 



-o'(3jc2+1) *aSx-£ 

1-x^ 

^ g(-2.y-tt) (3.x^-H) 



( 2 . 12 ) 



Dividing both sides by 3x*+l and simplifying yields 
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X[3(X+3)x*+X-l] _ _/ . «6x 

3}?TT 

(2.13) 

^(3x^+1) ^ a (-2x-o) 

1-^ l-x^ 



Next, we consider components of (2.13) that have 3x*+l in 
the denominator. The idea is to select a so that 



X [3 (X+3) jf^+A-1] -a6x ■ 0 mod(3jf*+l) 



This condition on a is essential in removing the singularities 
at ±(-3)‘*. For now, set a and B so that 



a=px , S=g 



where p and q are constants. Substituting px, p, and q into 
(2.12) for Q, a' and B, respectively, and simplifying results 
in 

(-3X^-9X)x^+ (2X*+10X)x^+ (X*-X) « 

{-3p^-9p-9b) x*+ (-p^+2p-6b) x^+ (-p-b) 

Next, we equate the coefficients. Thus, p = 2X and q = 
-X(X+1) which implies that a = 2Xx and B = -X(X+1). 
Substituting back into (2.11), it follows that Y = (X+l)(-2x). 
Therefore, the system (2.10) takes the form 
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(1-x^)*' = 2Xx^ + (1+3x^)T 

= -X(i+1)* + (A+l)(-2x)f 



-Kx<l 



A simple substitution of (A+l)f = f implies the following 
equivalent form 



Hence, (2.14) is a system of first order differential 
equations that is equivalent to (2.8). 

C. PROPER ANALYSIS 

The Priifer analysis is used to study the oscillatory 
behavior of trajectory solutions in the corresponding phase 
plane [Ref. 2:p. 312]. It is a very powerful method that is 
extremely useful . The objective of this section is to use 
this method to establish the theorems listed below for general 
systems of differential equations, and then apply these 
theorems to (2.14). These theorems for systems are adapted 
from the analogous theorems for Sturm-Liouvil le differential 
equations [Ref. 2:pp. 312-353]. 

Theorem 1. (STURM SEPARATION THEOREM) 

If (S f)^ and (i 7)2 are linearly independent solutions of 
the system 



(1-x^)#' = 2kx9 + (X+1) (l*3x^)W 
(l-x^)W' = + (X+1) (“2x)¥ 



-Kx<l . 



(2.14) 
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- (I f] ($) 



-KX<1 



and 

4B(-j?)-tY-al*>0 

where a, B, y, and Q are polynomials in x, then must vanish 
at one point between any two successive zeros of f 2 • other 

words, the zeros of ^2 

Theorem 2. (STORM COMPARISON THEOREM) 

Let (i f) and ($ i\) be nontrivial solutions of the systems 

( 1 -^) (1^) • (I ?) (I) “d ■ (I f) (n) 



where a, 3, y, Q. a, U, y, and Q are polynomials in x, and 

4jS(-C>) -[y-e]2>0 
0 <S< jg 4S ( -S) - [y -o ] *> 0 

4(jg-B) (-((?-]:» )-(y-Y-(a-o)]»>0 

Then, 5 vanishes at least once between any two zeros of i, 
unless (5 t)) is a constant multiple of (§ f). 

Theorem 3. (OSCILLATION THEOREM) 

Consider the associated Priifer differential equation 
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(l-x*)^=A.Scos*6+ [ (A.+1) Y-A.O] sinOcosO- (X+1) j?sin*0 
obtained from the system of differential equations 

- (il (I) 

where a, 0, y, and Q are polynomials in x. If 

4jg(-C))-[Y-a]2>o 

and B( -1) = X, such that -x/2 < x < x/2 for each X, then the 
solution 0 is a continuous and strictly decreasing function of 
X for fixed x on -Kx<l. Moreover 

lim 0 (x) =-«> lira 0 (x) =«> 

X-*O0 X-*—* 

for -Kx<l. 

Theorem 4. Consider the system of differential equations 
(2.15) which satisfies the Oscillation Theorem. If 

0(-l)“X, aJJ-o»<X<«> 

0(1)=O-J3H, i2=0, ±1, ±2, . . . 

A « 



then (2.15) has an infinite sequence of real eigenvalues 
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with 



lim 

J3-»o» 



and 



lim An=-«» 

j 3 -»-oo 



The eigenfunction belonging to the eigenvalue Xj has exactly 
n zeros in the interval -Kx<l, and is uniquely determined up 
to a constant factor. 

1. Phase Plane 

The objective of this section is to convert to polar 
coordinates the following general system of differential 
equations 



(1-jc^)*' = o® + 

“ jS® + 



(2.16) 



where a, B, y , and Q are polynomials in x. By introducing 
polar coordinates 

♦=r(x)cos0(x) , T=r (x) sin0(x) 

where 6 is positive in the counterclockwise direction, the new 
dependent variables become r and 6. 

Expressing r and 6 explicitly as functions of S and f 

leads to 
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0«airctein^ 



and 



(2.17) 



jr2^^2+Jf2 ^ (2.18) 

Then, by differentiating both sides of (2.17) and 
(2.18), an equivalent system of differential equations for 
(2.16) is presented as functions of r and 6 . Differentiating 
both sides of (2.17) yields 

. (2.19) 



Since 



sec*0 - l+teU 1^0 




^ 2+^2 

“<F" 



then (2.19) becomes 



d0 _ 
~dx " 



( 2 . 20 ) 



By substituting from the general system (2.16), 
equation (2.20) takes the form 

dx ^2 
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Substituting in polar coordinates, the differential 
equation for 6 becomes 

(l-x^)-^=ficos^0+ (y-a) sin0cos0-j?sin^0 (2*21) 

Finally, we differentiate both sides of (2.18) and 
follow similar steps. The differential equation for r then 
becomes 

(l-x*)^=r [aco8^0+ (j?+5) sin0cos0+ysin*0] . (2.22) 

Therefore, (2.21) and (2.22) represent an equivalent 
system to (2.16). 

The advantage of this system over (2.16) is that 

(2.21) is a first order differential equation of only 6. This 
allows the study of the oscillatory behavior of the solution 
trajectories independent of r. To find r for any solution 0, 
one solves (2.22) directly. 

We also note that r is always positive for nontrivial 
solutions. This can be seen by examining the solution of 

(2.22) [Ref. 2;p. 312]. 

Next, we recall the following result from elementary 
algebra. A function 

F{x,y) = ax*+jbAy+cy* 

is negative definite if and only if 
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a<0 and 4ac - jb^>0 



and positive definite if and only if 

a>0 and 4ac - jt>*>0 



Next, we apply this result to (2.21). If we can show 

that 



4jg(-£))-[Y-o]2>0 



then equation (2.21) is either negative definite or positive 
definite, depending on the sign of fl. If equation (2.21) is 
positive definite, then the change in the angle 6 is always 
greater than zero. This means that the angle 0 is always 
increasing. Since our polar coordinate system is set up so 
that the positive direction of rotation is counterclockwise, 
then the solution trajectory of 6 always rotates 
counterclockwise. If equation (2.21) is negative definite, 
then the solution trajectory of d always rotates clockwise. 
In addition, since the change in the angle 8 is either 
positive definite or negative definite, then the direction of 
the solution trajectory 8 can never change. Therefore, we can 
conclude that once a solution trajectory crosses an axis it 
can not cross back. 

Also, since i = rcos8, then the zeros of $ occur when 
8 = x/2+nx where n = 0, ±1, ±2, .... 



17 



Consider the system of differential equations (2.14). 
The associated Priifer equation is negative definite for X>0 , 
and positive definite for X<0. However, for (2.14), it is 
easy to determine the direction of rotation across the #-axis 
and ¥-axis without the Priifer equation. First, we examine the 
case when X>0 by directly substituting into (2.14). 

By setting $ = 0 we can determine the direction of 
rotation along the ¥-axis. Also, by setting f = 0 we can 
determine the direction of rotation along the §-axis. Hence, 
if we substitute these values back into (2.14), then the 
corresponding phase plane for X>0 takes the form in Figure 1. 




Figure 1. Phase Plane Analysis Along f-axis and §-axis 
for >l>0 
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clearly, the solution trajectories can cross the f- 
axis and f-axis only in a clockwise direction. In the case 
when X<0, the trajectories rotate in the counterclockwise 
direction . 

2. Sturm Separation Theorem 

Suppose ($ and ({ are two linearly independent 
solutions of the system 

(!') “ (I ?) (?) 



and 

4jg(-C>)-[y-«]*>0 

where a, B, y , and Q are polynomials in x. 

There are two cases to consider, one being when B<0 
and the other when B>0. First, we consider the case when B<0 . 

We examine the area "A" of the triangle created by 
joining the points (5 ?)^, (§ 1)^, and the origin in the (*,T) 
plane. Figure 2. 

By inspection, the area "A" is equal to zero only when 
(i T)^ and (§ T)2 are on the same straight line through the 
origin. This can occur with (I f)j and (5 f >2 both in the same 
quadrant, or in opposite quadrants. However, this can never 
occur since (§ f)^ and (5 f )2 are linearly independent, so "A" 
has positive area. 
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¥ 

/N 




V 



Figure 2. Area Principle Analysis 

Recall that two nontrivial solutions of (2.23) are 
linearly dependent if and only if they are located on the same 
straight line through the origin in the (*,f) plane. Hence, 
the angle, tt, between the linearly independent solutions, 

(i and (5 ¥> 2 / is always 0<«<5t. In other words, (S T)j and 
(f f )2 are always in the same rotating half plane and never 
cross each other. 

Analytically, we can show that the area "A” is never 
zero by the following method. Using the area principle, we 
calculate the area "A" using the following formula 



20 






= % 



(2.24) 



1 Y, 

1 #2 ^2 
10 0 



*2 ^2 



The change in the area is 

^ - V4 [♦',», . - *,¥11 

Substituting in f'j, ^1' ^2 from (2.16) and 

simplifying gives 

^-\4(«n) (*iTa - *2»i) 



Substituting into this equation from (2.24) results in 



dx 



(«+y)-A 



Finally, by solving for "A", the following formula is derived 
to determine the area "A" at any point -l<x<l; 

'(o+Y)dt| 



A“ic exp 



(A 



where k = A(-l) is the initial area. 

Since a and y ^i^e polynomials and -l<x<l, the area is 
finite. Moreover, A(-l) can not equal zero since the two 
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solutions are linearly independent. Hence, the area "A" can 
never be zero. 



Next, we recall that i3<0 and 
4jS(-C>)-[y-o]2>0 

Therefore, d6/dx is negative definite which implies that all 
solutions rotate in a clockwise direction. Also, once a 
solution crosses the f-axis, it can never cross back. 
Finally, since the area can not equal zero, the angle o 
between (f f)^ and (§ f )2 is always satisfies 0<o<x. 

Therefore, as the solutions rotate clockwise in the (i,f) 
plane, (f T)j and (f ¥>2 always cross the ¥-axis in the same 
order. Therefore, the zeros of and §2 interlace. The 

analysis for the case when B>0 proceeds in exactly the same 
way except that the solutions rotate counterclockwise, q.e.d. 
3. Sturm Comparison Theorem 

Let (§ ¥) and ({ n) be nontrivial solutions of the 

systems 

■ (fi ?) (») 

(2.25) 

“ (If) (n) 



respectively, where a, B, f, Q, a, TJ, y, and Q are polynomials 
in X, and 
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4£{-C>)-[Y-a]*>0 
o<B<£ 4S(-5) - [Y-«] ^>0 

4 (B--B) (- (C>-2>) ) - [Y-Y- («-«) 1 ^>0 

Again, there are two main cases to consider. One case is when 
B-ir > 0 and the other when B-'E’ < 0. The second case is the 
same as the first except that the solutions rotate in the 
counterclockwise direction. Hence, it suffices to consider 
only the case when B-E* > 0 and that all solutions rotate 
clockwise. The objective is to analyze the two solutions in 
the (J,f) plane through the use of the Priifer equation. 

The associated Priifer equations for (2.25) are listed 
below, respectively. 

HR 

-^=&:os^0i+ (Y~«) sin0^cos0j^ - (?sin^0j (2.26) 

(l-x^) « Bcos^ 02+ (Y“^) sin02cos02 - ^in^02 (2.27) 

By subtracting (2.27) from (2.26), the resulting 
equation indicates the rate of change of the angle formed 
between (S f) and (( i|). In addition, by superimposing one 
solution on top of the other, with 8 ^ = 82 (Figure 3), then the 
equation becomes 
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(1-x^) (0i"02) * (jS-S) cos*0i+ [Y“Y“ (ot-o) ] sln0^cos6^ 

- iO-D) sin^i 

The resulting equation is a Prvifer equation. This 
equation facilitates the analysis of the solution vectors 
({ f) and (f t)} by examining the behavior of the two angles, 
and 02 . 

Since fl-B" > 0, and 

4(S-£) (-(|?-0))-[Y-Y-(a-«)]^>0 

then, d/dx(0j-02) > 0. This means the rate of change of the 
angle between (f f) and t|) is always increasing, and (( t)) 
"runs away” from (f f ) , Figure 4. 

Note that (5 f) and n) are solutions of two 

different systems of differential equations. The two solution 
may be linearly independent or dependent. Hence, one solution 
could "lap" the other solution one or more times. 
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T 




Figure 3. Superimposition of Solutions 

T 




Figure 4. One Solution "Runs Away” from the Other 
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Recall that the zeros of i and ( occur when (i f) and 
({ q) respectively cross the f-axis. Therefore, as the 
solutions continue to rotate, between any two zeros of i there 
is at least one zero of Note that, if (S T) = c(( i)), for 

all -l<x<l, then the zeros occur simultaneously. Hence, ^ 
vanishes at least once between any two zeros of i, unless 
i|) is a constant multiple of ({ T) . q.e.d. 

4. Oscillation Theorem 

We consider the system of differential equations of 
the form 







Xa (X+l)ff 
XS (X+1 )y 




-l<x<l 



(2.28) 



where a, IB, y, and Q are polynomials in x, and 

4£(-C>)-[y-e]»>0 



The system (2.28) is basically the same as (2.23) 
except now we have introduced an eigenvalue, X. The objective 
now is to vary X thereby to study the zeros of f. The Priifer 
equation for (2.28) is 



(l-x*) -^i=)ljScos*0i+ [ (X+1) y-JLo] slnO^cosB^ 



- (X+1) 



(2.29) 
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As before, the zeros occur only on the T-axis which is 
the same as saying the zeros occur at 6 = x/2+nx for n = 0, 
11, 12,.... Next, we impose endpoint conditions on 0 such 
that 



0(-l)-T, -^<X<^ 



for all-«><X<«> 



0(l)=a-ii*, n=0, ±1, ±2, . . . 

^ A 



(2.30) 



so that the angle t is as close to x/2 as possible. Note that 
the endpoint conditions are computed from tan6 = f/5. 

Now, we consider a second system 

- (m inti!?) (^) 



where u is an eignevalue. The Priifer equation for this system 
is 

dB 

(l-x*) -3 -^=iiScos 202+ [ (|i+l) Y“ 1 »o 3 sln0,cos02 

dx 2 »- I ^2.31) 

-(ji+l)Osin*02 

Suppose that x is fixed and the differential equation (2.31) 
satisfies the following endpoint condition 

0(-l)»t, for 

2 2 
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Following the same basic proof as in the Sturm Comparison 
Theorem, we subtract (2.29) from (2.31), and equate to 02 
at 0j(-l) = 02(-l) = T. Hence, 



(j»-A) [Scos^Q^* (y-a) sin0iCOS0i-(?Bin*6i) 



(2.32) 



Since 

4j6(-C>) -[y-«)^>0 

then the discriminant of (2.32) is 

(»i-A)M4ig(-C>)-(Y-a)^J>0 . 

Therefore, d/dx(©2-02) is either positive definite or 
negative definite. Accordingly, there are several cases to 
consider. However, by examining one case carefully, all other 
cases follow directly, with only the direction of rotation 
varying. Thus, it suffices to show only one case. 

Let X be fixed. Suppose that 0<p<A which implies that 
U-A < 0. Now, consider the case when fl<0. 

Since fl<0 and p-A<0, then d/dx(02-0j) > 0. Hence, by 
considering the conditions that 0<)i and 0<A, it is implied 
that all solutions rotate in the clockwise direction, and 
(# f) "runs away" from (5 t)) as A goes to positive infinity. 
Figure 5. 
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▼ 




Figure 5. (# T) "Runs Away" from t|) , A,>0 

If A,<0, then u~A > 0 which implies that 
d/dx( 62 “®i ) < 0- Therefore, when A.<0, all solutions rotate in 
the counterclockwise direction and (5 ?) "runs away" from ($ 
f|) as X goes to negative infinity. Figure 6 . 

In summary, as A. goes from -« to », the angle 6 ^ goes 
from « to -« in a strictly decreasing manner. Hence, for 
fixed X, 6 is a strictly decreasing function of the variable 
X. In addition, 

lim 0“-o9 , liin 

A— ^ A— <• 



q.e.d. 
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¥ 




Figure 6. (i ¥) "Runs Away" from (5 H) / ^<0 

5. Sequence of Eigenvalues and Interlacing Zeros 

Consider the system of differential equations (2.15) 
which satisfies the Oscillation Theorem. Now we continue with 
our analysis and maintain the conditions of the previous 
proof. Namely, let 0<u<X and 

^oz all 

A A 

6(l)*o-rm, O‘‘0,±l,±2, . . . 

A A 

where x is as close to x/2 as possible. 
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Now, we let y be fixed and B<0. Note that a value for 
X which satisfies both endpoint conditions is an eigenvalue 
for the system of differential equations (2.28). Let Oj(-l) 
= X. Then, as X increases, the solution trajectory of 
rotates in a clockwise direction. Since d6j/dx < 0, the 

trajectory of eventually crosses ©^(l) = o-nx (n = 0). 

Hence, the first time that the trajectory of Oj crosses the 
line 0j(l) = o-nx is when n = 0. We call this eigenvalue Xg . 
Therefore, when X = Xq , the solution trajectory of ©j does not 
cross the T-axis. Hence, there are no interior zeros for the 
corresponding eigenfunction ig of Xg , Figure 7. 




Figure 7. Interior Zeros of fg 
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Recall that for 



each Xj , the corresponding 
eigenfunction is determined by the formula 

O^-rcosO 



Now, as X continues to increase, the solution 
trajectory of continues to rotate in the clockwise 

direction. Eventually, the trajectory of again crosses the 
line 6^(1) = o-nx when n = 1, or 6^(1) = a-%. We call this 
eigenvalue Xp Therefore, when X = X^, the solution trajectory 
of crosses the f-axis one time. Therefore, the 

eigenfunction of X^ has one interior zero. This process 
continues as X increases to «. Similarly, for X<0, the 
solution trajectory of 6^ rotates counterclockwise and this 
process continues as X decreases to -®. 

In svunmary, as X goes to ±«, each time that the 
solution trajectory of 6^ crosses the line 6^(1) = o+nx, the 
corresponding value for X becomes an eigenvalue of the system 
(2.28). Hence, there is an infinite sequence of Xj, for n = 0, 
±1, ±2,... for which the second endpoint condition is 

satisfied. Therefore, there are eigenvalues, Xj, such that 
. . . <X.j^<Xj<X 2 < . . . , and the associated eigenfunction belonging 
to Xj has exactly n interior zeros in the interval -Kx<l. 
Moreover, it follows that 



llmX„ 

Jl-*oo 



and 



limX^*-«o 

n-*-» 
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In addition, any two solutions which satisfy the same 
initial conditions are linearly dependent. Therefore, each 
eigenvalue uniquely determines an eigenfunction up to a 
constant factor. This completes the proof of Theorem 5. 

For the system of differential equations (2.14), the 
discriminant of the associated Priifer equation is positive. 
Therefore, all of the previous theorems apply. 



D. ORTHOGONALITY AND COMPLETENESS 

The objective of this section is to show that in a Hilbert 
space the eigenvectors of (2.14) are complete and orthogonal 
with respect to a weight function. [Ref. 3:pp. 344-353] 

1. Orthogonality. 

The system of differential equations (2.14) is 
equivalent to the following form 

U-x“)($;) - (♦) * (° 1!|^) (I). (2.34) 

For convenience, represent (2.34) with the notation 
(l-x*)Zi-XAZi+BZj^ (2.35) 

where Zj is the eigenvector (5 f) associated with the 
eigenvalue X. A and B are the 2x2 matrices given below 
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A 




2x l+3x* 
-1 -2x 



and B - 



0 i+3jc^ 
0 -2x 



) 



In addition, we let L stand for the linear 



differential operator, such that, when L operates on Z then LZ 
is defined by 



The general strategy is to find the adjoint operator. 



we integrate both sides to establish orthogonality. Finally, 
matrix algebra can be used to find the corresponding weight 
function . 



W^(I2) - (l-x^)W^'-XW^AZ-ff^Z 

W^]'-[(l-x^)hr'n'Z-XW^AZ-W^Z 

Then, by requiring the exact derivative 

lf^{UZ)-(L*W^Z‘-^Ul-x^)W^] , 

dx 

it follows that 

Z=-{ [ (1-x^) ' Z+XRT’AZ+J^^BZ} 



LZ= il-x^) Z'-XAZ-BZ=0 



L|, such that W^(LZ) - (L*W^)Z is an exact derivative. Then, 



Let W be a 2 X 1 column vector such that 
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Hence , 



L*l>^= -U ( 1 -x^) W '+ 

For general u, the adjoint differential equation is 

or (2.37) 

[ (1-x*) W] ' - 



Next, we return to (2.36). Substituting w'(LZ) and 
(lV)Z into the equation and simplifying yields 



W^(UZ)-(L*W^Z‘^[W^(l-x^)Z] 

ax 



which simplifies to 



^’•(l-x*) Z4 ] > Z--^ IW^(l-x^) Z] 



Substituting in (2.35) and (2.37) results in 



(A-|i) W^AZ‘-^ [W^(l-x^) Z] 
dx 



Now, we integrate both sides. This leads to the following 
orthogonality relationship 



{X-H)pW^AZ^(W^{l-x^)Z)\lj, »0 



(2.38) 
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If the integrand is well defined, then the 
eigenvectors W are orthogonal to those of Z, with respect to 
the matrix weight function A. Next, we determine the nature 
of W. 

Expanding (2.37) leads to 

(l-x^) (2.39) 



We write W in the form H = JZ where Z is a 2 x 1 column 
vector, and J is a 2 x 2 matrix of the form 



J- 




Then, (2.39) becomes 



( 1 -X* ) JE'=2xJE 



Then we, premultiply both sides by J*^. Note that 
J*^ = -J. Thus, 

(l-x*)E'=2xE+»iJA^JS+JB’'JE . (2.40) 



Since 



J) (?? 



3x^+1 

-2X 




A 



and 
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(2x 3jf2+l\ 

i 0 0 j 



then (2.40) simplifies to 

(l-x2)E'=^iAE + (B+4jd)E (2.41) 



where I is the identity matrix. Next, letting 

E-(l-x2)-2Q and 

E-4x(l-x*) "^Q+ (1-x^) ‘^Q 

where 0 is a 2 x 1 column vector and substituting into (2.41) 
gives 

(l-x2)Q'=|iAQ+BQ . (2.42) 

So, Q = Zp where Zp is the eigenvector corresponding to the 
eigenvalue \x. Back substituting results in 

JE » t7(l-x*)-*Q - JH-x^)-^Z^ 

If 1 # y, then the orthogonality relationship (2.38) 

becomes 

Since A ^ y, then X-y * 0 which implies 
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When X = u, then 




J'^k 



Zx dbc ^ 0 



(2.43) 






since the integrand is always positive. In addition, since 
the integrand is a polynomial, then it is finite. This 
implies Zj is square-integrabl e with respect to the matrix 
weight function 



in (2.43), and all the remaining entries are polynomials. 

Next, we determine the sequence of eigenvectors. 
Clearly, X = 1, 2, 3,... all yield solutions. This is based 
on the definition of X such that X = n-1 where n = 2, 3, 
4,.... Now, we consider X = -1, -2, -3.... Since 



for the Legendre polynomials, the following relationships 
hold. 



J-*A 






Initially, the (1-x*)*^ term may seem to imply that the 
integrand is not well defined at ±1. However, (1-x*) cancels 
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^X*l “ “ ^-X-7 



( 2 . 44 ) 



*X*2 “ -P-Ca^a)*!) ^-X-3 



(2.45) 



Substituting (2.44) and (2.45) into the original 
definition of j yields 



♦n.l = 



P, P, 






« 



n,l 



Pi P2 



P-X-2 P-X-1 



for X 2 1. Equivalently, for X s -1, 



Kx 



Pi P2 

Px-2 Px~2 



This form provides the solution to the original 
differential equation (2.8). Therefore, there exists a 
sequence of eigenvectors 








where (f f) has polynomial entries of degree < m+S. 

2. Completeness 

Finally, we determine that the sequence of continuous 
eigenvectors (f f)j for m = 0, ±1, ±2, ±3,..., which are 
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orthogonal and square-integrable with respect to a matrix 
weight function, is complete. That is, we show that every 
square-integrable function F = (1-x*) (f g)* can be expanded 
into an infinite series such that 

'‘•'’(J) ■ ,1 “• («). ■ 

Recall that there is a factor (1-x*) which cancels from both 
sides of the equation. 

Consider the coliunn vector (p q)^ where p and q are 
polynomials of degree < m+3. There exist constants c.^ and c^ 
such that 

(5) -I I 



is a column vector with polynomial entries of degree <m+2. 
This follows by equating the coefficients of and solving 
the resulting system of equations for the coefficients c.^ and 
c^. Since (J ¥).^ and (f T)^ are linearly independent, then c.^ 
and Cg exist. By mathematical induction, one can find 
constants C|^ for k = -m, ... , -1, 0, 1, ... , m such that 
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Finally, we let P = (1-x*) (f g)* where f and g are 
any continuous functions on a finite closed interval -l<x<l, 
and let Ep €2 > 0 . Then, by the Weierstrass Approximation 
Theorem, there exist polynomials such that 

and 



for all X on [-1,1]. Hence, 



(^) - («) 



( 



(l-x^) t-P 
(l-x^)g-q 




^ max 




which implies 



/ X 


A ® /A\ 




) ■ '• (?) 



i max (€^,€2} 




Therefore, F is approximated uniformly, within an 
arbitrarily small distance 6 > 0, by a linear combination of 
the (f f ) . Hence, the sequence for m = ... , -1, 0, 
1, ... is a uniformly convergent sequence, which implies the 
sequence is mean square convergent on -l<x<l with respect to 
an integrable positive weight function. Recall that (1-x*)^ 
cancels in the integrand, leaving a positive definite matrix. 
Therefore, the sequence is complete. [Ref. 2:p. 353] 
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In addition, the constants C|^ are determined in the 
following way 



r Zj— dx 
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Ill . 



RESULTS AND CONCLUSIONS 



A. RESULTS 

We have shown in this thesis that there is a special 
function which takes the form 






(n>k) k=0,l,2, . . . 

(-Kx<l) 13=1,2,3,... , 



(3.1) 



where the elements of the determinant are Legendre 
polynomials, and possesses the essential Sturm-Liouvi 1 1 e 
properties. In addition, several interesting results are 
obtained and listed below. 

1. Hypergeometric Functions 

The Legendre differential equation is a special case 
of the hypergeometric differential equation. This suggests 
that a determinant of the form (3.1), with other 
hypergeometric functions as elements, might possess similar 
properties . 

Recall that the recursion formula 

Pn*l - ^(n^) P'n M 

and the Legendre differential equation are used extensively in 
the initial development of (3.1). In fact, it can be shown 
that since the hypergeometric function can be expressed in a 
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recursion formula of the form 



F(r+1 ; x) = A(r i x)F(r ; x) + B(r ; x)F'(r ; x) , 

then the determinant (3.1), with other hypergeometric 
functions as elements and appropriate boundary conditions, 
possesses the same Sturm-Liouvi 1 1 e properties. 

2. Turan's Inequality 

The proof of Turan's inequality is a straight forward 
application of the general case (3.1). Setting n = k+1 and 
substituting back into (3.1) results in a determinant that is 
of Turan's form. Recall, from the previous analysis of (3.1), 
that there exists a X such that (3.1) has no zeros in the 
interval -Kx<l. Hence,. Turan's inequality is established. 

3. General Systems of Differential Equations 

In general, boundary value problems of systems of 
differential equations are not well-posed. Until now, very 
little has been established concerning systems with polynomial 
coefficients. In Chapter II, several theorems are developed 
and proved for such systems. Hence, if the conditions of the 
theorems are satisfied, then the system is said to be of 
Sturm-Liouvi 1 1 e form. 
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B. CONCLUSION 



The analysis of this one problem presents several possible 
avenues of further research. Only the first step was taken in 
a direction which could result in establishing several new and 
substantial results. For example, the next step would be to 
establish the case when (3.1) has other hypergeometric 
functions as elements. Then, one could continue the research 
on Sturm-Liouvil 1 e systems of differential equations by 
examining the possibility of expanding the necessary and 
sufficient conditions for such problems. Research on this 
subject is ongoing. 
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APPENDIX 



I. PROOF FOR GENERAL CASE - *n|j(x) 

A . PURPOSE 

In this appendix we briefly outline the procedure for 
establishing the Sturm-Liouvi 1 1 e properties for the general 
case, The basic procedure is very similar to that 

illustrated in Chapter II. Thus, only the important steps are 
highl ighted . 



B. DEVELOP DIFFERENTIAL EQUATION 

By following the procedure in Chapter II, it is easy to 
establish the following relations for $jjjj(x), $'u|j(x), and 



®=- (1-x^) 



n+l 



pA' 

k+1 



( 1 ) 



^/- (n-k) (1-x^) 
■ (71+1) ik+1) 



p'A 



(n-k) P^Pt, 



( 2 ) 



<!>"= (n-k) 



(n-k+1) 

n+l 






2x 

(71+1) (k+1) 



pA* 



(^~'^'*’l ) p p! 

ic+1 “ * 



(3) 



Equations (1), (2), and (3) lead to the differential 
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equation 









(1-^) 

(ic+D* 



(Pi) 2 






2 

ic+1 






2x 

(Ic+D* 




*4 



(4) 



(fl-lc) 



(fl-ic+l) Pjt*+ 



2x „ (iC-J3+l) /r/\* 



♦ =0 



where the general solution is 



C\ 



P,. 

Pa 



Pt.i 






Pi. 

Oa 



Pt*i 

Oa*i 



Applying appropriately translated boundary conditions 
eliminates all logarithms, , in the general solution. Now, 
setting X = n-k-1 and letting 



(?'= 



2^ 

Jc+1 




( 6 ) 



then (4) takes the form 



(1-x*) (1-.X*) -^*4 (X+1) 

0 



x+. 



2P*(P, 



* lc+1 



Pi) 



*-0 



(7) 
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C. REDUCTION TO A SYSTEM OP FIRST ORDER DIFFERENTIAL 
EQUATIONS 

Next, we reduce (7) to a system of first order 
differential equations of the form 



“ «4> + -Kx<l 



( 8 ) 



where a, B, y, and Q are polynomials in x. 

First, note that Q and Q' have no zeros in common. This 
implies that all polynomial coefficients of J, , and are 
reduced to lowest terms and only simple zeros occur in the 
denominator. The proof is as follows. 

If Q' =0 , then either P/ =0 or 

* Pt~0 . (9) 



If P]j'=0, then Q=0 only when P|[ = 0. However, it is easy to show 
that P]j and Pjj' are relatively prime. Hence, when P|j'=0, then 
Q * 0. 

If (9) holds, then by solving for Pj^ and substituting back 
into Q it follows that Q=0 only when Pjj'=0. But, by the first 
case, this can not occur. Therefore, Q and Q' have no zeros 
in common. 

By transforming (8) into a second order differential 
equation and equating coefficients of f* and i' it is easy to 
establish that 
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Y = -2x-a 



( 10 ) 



Equating the coefficients of 5 and substituting (10) for y, 
results in the following equation. 



QS _ tt(-2x-tt) _ . 

1-A* 1-JC* Q 



-A(A+1) 



X+1 

0 






( 11 ) 



Now, we isolate and examine the components of the equation 
that have Q in the denominator. The idea is to select a such 
that 



aO'- (A+1) 






= 0 mod(C>) 



Substituting (6) for Q' implies 



k+1 



l^(;^l^*Pt)]-(^+l> ] =0mod(£)) 



However , 



X 

k+1 



Pi * Pt 



( 12 ) 



and Q are relatively prime. This is shown by the following 
argvunent . 
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Set (12) equal to zero and then solve for P^. By 
substituting back into Q it follows that (12) and Q have zeros 
in common when P]j and P|j' have zeros in common. However, Pjj 
and P|j' are relatively prime. Hence, (12) and Q are 
relatively prime. 

Therefore, since (12) and Q are relatively prime, then 
there exists a polynomial M such that 






(13) 



By examining (13) it follows that 

(X+1) = 0 mod(Pi^) 

Moreover, since Pjj and P|j' are relatively prime, then there 
exists a polynomial N such that 

(X+l) + MPj^ - . (14) 

However, by the Euclidean Algorithm for Polynomials, there 
exist unique polynomials and such that 

^Pj^ + //iPi “ 1 . 

From (14), this implies that 

Af - -(X+l)Af, and N= (X+DA^i 



50 



Hence, there are unique polynomials M and N that determine 
a such that all the zeros of Q divide out. 

Finally, choose B to satisfy (11). 

D. STURM-LIOUVILLE PROPERTIES 

Although it is somewhat tedious, it can be shown that 

4jS(-i?)-[Y-a]2>0 

Therefore, all of the theorems previously established 
apply. In addition, the properties of orthogonality and 
completeness follow in the same manner. 
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